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(57) The invention relates to an apparatus and a 
method for determining an approximate value for the 
stroke volume and the cardiac output of a person's 
heart. The apparatus and method employ a measured 
electrical impedance, or admittance, of a part of a per- 
son's body, namely, the thorax. This part of a person's 
body is chosen because its electrical impedance, or ad- 
mittance, changes with time as a consequence of the 
periodic beating of the heart. Accordingly, the measured 
electrical admittance or impedance can provide infor- 
mation about the performance of the heart as a pump. 
The apparatus and method are not restricted to the tho- 
rax, but it is also possible to insert electrodes located on 
a catheter/probe into the esophagus of a patient. 
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Description 



[0001] The invention relates to an apparatus and a method for determining an approximate value for the stroke 
volume and the cardiac output of a person's heart. The apparatus and method employ a measured electrical impedance, 
5 or admittance, of a part of a person's body, namely, the thorax. This part of a person's body is chosen because its 
electrical impedance, or admittance, changes with time as a consequence of the periodic beating of the heart. Accord- 
ingly, the measured electrical impedance or admittance can provide information about the performance of the heart as 
a pump. 

10 RELATED PRIOR ART 

[0002] In 1966, Kubiceketal. were the first to design a clinically applicable device, capable of determining the stroke 
volume (SV) by non-invasive, electrical means. The Kubicek method is disclosed in the article by Kubicek et al., De- 
velopment and Evaluation of an Impedance Cardiac Output System, Aerospace Medicine 1966, pp 1208-1212, and 

15 in US patent 3,340,867 which are both incorporated herein by reference, (see also US patents 5,178,154 to Ackmann 
et al., 5,316,004 to Chesney et al, 4,953,556 to Evans, 5,685,316 to Schookin et al, 5,505,209 to Reining, 5,529,072 
to Sramek, 5 ; 503,157 to Sramek, 5,469,859 to Tsoglin et al., 5,423,326 to Wang et al, and 5,309,917 to Wang et al.) 
[0003] When a tetrapolar array of circumferential band electrodes is placed at the base of the neck and about the 
circumference of the lower chest, at the level of the xiphoid process, and a constant magnitude alternating current (AC) 

20 is injected through the upper cervical and lower thoracic band electrodes, a voltage, proportional to the thoracic elec- 
trical impedance (or reciprocally proportional to the admittance), is measured between the inner cervical and thoracic 
band electrodes. The portion of the cardiac synchronous impedance change, AZ(t), temporally concordant with stroke 
volume, was ascribed solely and uniquely to volume (plethysmographic) changes of the aorta during expansion and 
contraction over the heart cycle. 

25 [0004] In the article by Woltjer H.H. et al. (The technique of impedance cardiography. Eur Heart J 1977; 18: 
1396-1403), the Kubicek model is explained as follows. The aorta is considered a cylinder of length L, equal to the 
distance between the voltage sensing electrodes. The thorax, exclusive of the aorta, is considered a cylinder of length 
L, equal to aortic length, and of cross-sectional area (CSA), equal to the cross-sectional area of the thorax measured 
at the xiphoid level. The blood-filled aorta is assumed to have a constant specific electrical resistance equal to that of 

30 stationary blood, p. The thoracic encompassing cylinder is assumed to be homogeneously perfused with blood of 
specific resistance p. The aorta and the thoracic encompassing cylinder are assumed to be analogous to parallel 
electrical conductors. 

[0005] It was accepted by Kubicek that, according to Nyboer (Nyboer J. Electrical impedance plethysmography. A 
physical and physiologic approach to peripheral vascular study. Circulation 1950; 2: 811-821), the portion of AZ(t). 

35 temporally concordant with SV. represented simultaneous inflow and outflow of blood over the systolic portion of the 
heart cycle. Thus, determining the area underneath the systolic portion of AZ(t) was assumed not to represent net 
volume inflow across the aortic segment under electrical interrogation. Thus, an extrapolation procedure was proposed, 
utilizing the maximum forward systolic slope of AZ(t). In order to compensate for aortic outflow, the maximum forward 
slope, analogous to peak flow, was stipulated to be constant throughout the systolic ejection interval. 

40 [0006] The maximum forward systolic upslope represents the peak, or maximum rate of change of impedance, i.e. 



[0007] Instead of measuring the slope directly, as proposed by Nyboer, Kubicek electronically differentiated AZ(t) 
into dZ(t)/dt. 

Thus, the peak systolic magnitude of dZ(t)/dt is 



[0008] In order to derive stroke volume (SV), Kubicek multiplied the peak rate of change of impedance by systolic 
flow time of the left ventricle, T LVE . 
[0009] According to Kubicek 
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sv = v 



I dt ) MAX 



iff 



T LVE =P- 



dZjt) 
dt 



MAX 



*LVK 



wherein Z 0 is the quasi-static portion of the measured impedance Z, and 
wherein 



< dZ{tf \ 
k dt j 



MAX 



is the peak value of the (inverted) first time-derivative of AZ(t), which corresponds to the maximum forward systolic 
upslope of AZ(t). Note that in this context, by peak magnitude, the maximum absolute amplitude is stipulated. In fact, 
during systole, the impedance decreases such that the sign of AZ(t) is negative. Hence, correctly stated, 



dZ(t) 

V dt ) 



MAS 



is the minimum of the time-derivative of AZ(t), i.e. 



dZ(t)) 



dt 



Jmin 



[0010] Furthermore, in the above formula, T LVE is the left ventricular ejection time, i.e. the time between opening and 
closure of the aortic valve, also referred to as systolic flow time. The volume 



V EFF =p — 



is the volume of electrically participating thoracic tissue (VEPT), wherein p is the specific resistance of stationary blood, 
which Kubicek assumed to be 150 £2cm, and L is the distance between the voltage-sensing electrodes which are 
applied to the neck and thorax. 

[0011] By virtue of rigid theoretical constraints, the Kubicek method, and its derivatives, consider volume changes 
in the aorta, i.e. plethysmographic changes, to be the sole contributor to 



dZjt) 

dt 



MAX 



[0012] Consequently, AZ(t) is assumed to represent the time-variable volumetric expansion and recoil of the aorta. 
Thus, its time-derivative, dZ(t)/dt, represents an ohmic equivalent of the rate of change of aortic volume. This would 
also imply that 



dZjt) 
dt 



MAX 
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measured in [QJs], is directly proportional to peak flow [mL/s] and peak velocity [cm/s], 

[0013] It is widely believed that the assumptions made in the Kubicek model are generally valid, i.e. that the increased 
aortic volume during mechanical systole leads to the decrease in the thoracic impedance. Since Kubicek assumed a 
directly proportional, i.e. linear, relationship between SV and 



dZ{t) 
dt 



J MAX 



times T LVE , it is usually believed that 



dZ{t) 
dt 



MAX 



is analogous and proportional to peak flow, or peak rate of change of aortic volume. Therefore, subsequent improve- 
ments focused only on a better definition and modeling of V EFF . 
[0014] For example, Sramek developed a formula according to which 

L 3 

Veff ~ 425 

(see US patent 4,450,527 which is incorporated herein by reference). 

[0015] In a later iteration, Sramek approximated L as 17 % of the person's height h. Thus, Sramek proposed the 
equation 



(o.nh) 



dZ{t) 



dt 



gy V v - A ' '*/ 9 \ _ /MAX t rp 



4.25 Z 0 



Bernstein (Bernstein D.P., A new stroke volume equation for thoracic electrical bioimpedance. Crit Care Med 1986; 
14: 904-909) introduced a factor 8 accounting for the person's weight deviation from ideal (as determined from the 
Metropolitan Life Insurance tables), corrected for blood volume, normalized to deviation from ideal body weight. Oth- 
erwise, Sramek's model remained unchanged, and Bernstein proposed the formula 



(dZ{t) 

sv = s (0-17*) 3 I dt /MAX T 



4.25 Z 0 



l.VK 



[0016] Despite these various efforts for improving the determination of the stroke volume, the stroke volume could 
not be correctly predicted across a wide range of subjects in health and disease states. 

[0017] In particular, in the following cases, the Sramek-Bernstein equation generally results in an overestimation of 
the true predicted stroke volume: 

children and healthy young adults; underweight individuals; tall, thin adults. 

[0018] According to Spiering et al (Comparison of impedance cardiography and dye dilution methods for measuring 
cardiac output. Heart 1998; 79: 437-441), the use of the Sramek-Bernstein equation generally results in an underes- 
timation of the true predicted stroke volume in the following cases: elderly adults, obese individuals; individuals with 
sepsis, acute lung injury or pulmonary edema; and during exercise. 

[0019] It is therefore an object of the invention as defined in the appended claims to provide apparatus and method 
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that determines the stroke volume accurately for individuals of all ages in health and disease states. 
[0020] The invention considers the absolute peak rate of change of impedance, 



( dZit) 



V dt J MIN ' 



to be the ohmic equivalent of peak aortic blood acceleration [ml_/s 2 ], or peak rate of change of aortic blood velocity. 
As a consequence, AZ(t), in earliest systole, is related to hemorheologic (blood flow) changes, not plethysmographic 
(volume) changes. Thus, the new apparatus can be described as an 'electrical velocimeter', or the method incorporated 
as 'electrical velocimetry'. 
[0021] Consequently, the measured value of 



must be integrated in order to obtain an ohmic equivalent for blood velocity. In summary, the invention mandates that 
the part of the previous art related to 



be changed. 

[0022] Hence, the apparatus and method according to the invention employ no underlying modeling or theoretical 
assumptions of the Kubicek, or any other subsequent, plethysmographic approaches. 

[0023] Furthermore, the approximation of SV according to this invention is not limited to the impedance method, but 
can be performed using the admittance approach. With 




cannot be implemented directly into SV calculation. Theoretically, 







and 



( drm _ i (dz{t)\ 



I dt Z 0 2 { dt ) MlN ' 



[0024] the SI/ is approximated according to 
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dY{t) 



sv = v, 



MAX 



[0025] With respect to the admittance approach further details are given towards the end of this description. 
[0026] It should be understood thatthe appended claims encompass the admittance approach (apparatus and meth- 
od) if the above-mentioned substitution of the impedance term is made. 

[0027] According to theory derived from basic science (and published as Sakamoto K, Kanai K. Electrical charac- 
teristics of flowing blood. IEEE Trans Biomed Eng 1 979; 26: 686-695; Visser KR. Electrical properties of flowing blood 
and impedance cardiography. Ann Biomed Eng 1 989; 1 7: 463-473; Lamberts R et al. Impedance cardiography. Assen, 
The Netherlands: Van Gorcum 1984; 84-85; and Matsuda Y et al. Assessment of left ventricular performance in man 
with impedance cardiography. Jap Circ J 1978; 42: 945-954), the change of blood resistivity, and the rate of change 
of blood resistivity, can be normalized for corrected flow time, FTC, 



where T LVE equals the left-ventricular ejection time (known also as systolic flow time), and is divided by the m th root 
of the cardiac cycle period T RR , where T RR equals the value for the RR interval (cycle time) in seconds and 0 < m < 1 .5. 
[0028] With V EFF defined as the effective volume of electrical participating thoracic tissue ([V EFF ]=m\), the stroke 
volume SV f according to the invention, is calculated from measured impedance according to the formula 



(or from measured admittance ) 
with 0.15 < n < 0.8 and 0 < m < 1 .5, 

and wherein C 1 is a constant which is necessary if n + m ^ 1 in order to adjust the units of the measured values in the 
formula such that the stroke volume is obtained in milliliters. C 1 need not have a numerical value different from 1 . 
[0029] A preferred case is that n = 1 - m. Then, 0-1=1. 
[0030] The most preferred case is n = m = 0.5. Then, 



[0031] Other objects, features and advantages of the invention will become apparent from the following description 
of a preferred embodiment of the invention. 

Fig. 1 is a schematic diagram depicting the apparatus, which uses on a subject's left and right side of the thorax, 
each, a tetrapolar surface electrode array, an alternating current (AC) source and a voltmeter. 

Fig. 2 is a schematic diagram depicting the apparatus of Fig. 1 , but limiting the measurement to the left side of the 
subject's thorax only. 

Fig. 3 is a schematic diagram depicting the apparatus of Fig. 2, using a tetrapolar surface electrode array across 
the subject's thorax. 
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Fig. 4 is a schematic diagram depicting the apparatus of Fig. 1 , but utilizing a common AC source for left and right 
side of the subject's thorax. 

Fig. 5 is a schematic diagram illustrating the application of an esophageal catheter / probe for determining SV by 
5 measurement of esophageal electrical bioimpedance (or bioadmittance). 

Fig. 6 is a diagram showing the variation of thoracic impedance with a ventilation cycle and with each heart beat. 

Fig. 7 is a diagram showing the variation of thoracic impedance with each heart beat when ventilation is sup- 
10 pressed. 

Fig. 8 is a diagram showing curves representing a surface electrocardiogram, the time-varying portion of the car- 
diogenic thoracic electrical impedance, AZ(f), and the time-derivative of this varying portion of the impedance, 
dZ(t)/dt. 

15 

Fig. 9 is a block diagram of the apparatus of the present invention with details of the different components of the 
apparatus; 

Fig. 10 is the block diagram of Fig. 9 in which the Digital Signal Processor is shown in less detail whereas the 
20 Microprocessor is shown in more detail. 

[0032] Principally, an alternating electrical field is applied to a thoracic volume forcing an alternating current (AC) to 
flow in parallel to the direction of aortic blood flow, i.e. vertically between neck and lower thorax. The current of known 
magnitude causes, in the direction of the electrical field, a voltage drop, which is proportional to the product of thoracic 

25 impedance and current applied. 

[0033] Fig. 1 schematically shows an apparatus according to the present invention, and its electrical interface with 
a subject 10. For the measurement of transthoracic electrical bioimpedance (or bioadmittance), a tetrapolar surface 
electrode array (with electrodes 12, 14 for AC application, and electrodes 16. 1 8 for voltage sensing) is applied to the 
subject's left side, and another tetrapolar electrode array (with electrodes 20, 22 for AC application, and electrodes 

30 24, 26 for voltage sensing) to the subject's right side. 

[0034] The left sided electrode array includes two current electrodes 12, 14, which are connected to an AC Source 
L 28, and two voltage sensing electrodes 16,18, which are connected to Voltmeter L 30. One voltage sensing electrode 
(16) is placed at the base of the neck, the other one (18) the lower thorax, at the level of the xiphoid process. The 
current electrodes (12, 14) are placed respectively, in the vertical direction, above and below the voltage sensing 

35 electrodes (16, 18). 

[0035] The right sided electrode array includes two current electrodes 20, 22, which are connected to an AC Source 
Ft 32, and two voltage sensing electrodes 24, 26, which are connected to Voltmeter R 34. One voltage sensing electrode 
(24) is placed at the base of the neck, the other one (26) the lower thorax, at the level of the xiphoid process. The 
current electrodes (20, 22) are placed respectively, in the vertical direction, above and below the voltage sensing 
40 electrodes (24, 26). 

[0036] AC Source L 28 and AC Source R 32 are voltage-controlled current sources (VCCS). Each VCCS provides 
an alternating current (AC), which is measured via Analog/Digital Converters 42, 44. Alternatively, the magnitude of 
the alternating current can be held constant, and the Analog/Digital Converters 42, 44 can be omitted. A Digital/Analog 
Converter (DAC) 36 provides an output that controls the AC Source L 28 and AC Source R 32. The Digital/Analog 
45 Converter (DAC) 36 itself is controlled by a Signal Synthesizer 38. The Signal Synthesizer 38 is implemented via a 
lookup-table in the memory of a Digital Signal Processor 40 as part of a Processing Unit 80 (indicated by dashed lines). 
Alternatively, a Direct Digital Synthesizer (DDS. not shown) can provide the functions of DAC 36 and Signal Synthesizer 
38. The Processing Unit 80 recognizes AC magnitude and phase of each VCCS. 

[0037] The voltages measured by the Voltmeters L 30 and R 34 do not only contain a signal caused by the AC 
so applied, but also a signal component from which an electrocardiogram (ECG) can be derived. The application of filters 
separates the AC related and ECG related signal components. The AC related signal component is proportional to the 
product of current applied and the impedance (which is unknown). In the case that the currents applied are of constant 
magnitude, the voltage V L 46 obtained by voltmeter L 30, and digitized by the analog/digital converter 48, is directly 
proportional to the unknown impedance of the left hemi-thorax, Z L (t) (or reciprocally proportional to the unknown ad- 
55 mittance, Y R (t)). With the AC magnitude held constant, the voltage V R 50 obtained by voltmeter R 34, and digitized by 
the analog/digital converter 52, is directly proportional to the unknown impedance of the right hemi-thorax, Z R (t) (or 
reciprocally proportional to the unknown admittance, Y R (t)). The Processing Unit 80 determines Z(t) by averaging Z L 
(t) and Z R (t), or Y(t) by averaging Y L (t) and Y R (t). Alternatively, the voltage V L 46 sensed between electrodes 16 and 
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18 (left side) and the voltage V R 60 sensed between the electrodes 24 and 26 (right side) can be summed, or averaged, 
prior to the voltmeters, requiring then only one analog/digital converter. 

[0038] A demodulation of the AC related signal component is required in order to extract the impedance related 
information from the AC carrier signal. Demodulation of the voltages obtained from the thorax is described, for example, 
by Osypka and Schafer (Impedance cardiography: Advancements in system design. Proceedings of the X. International 
Conference on Electrical Bio-Impedance (ICEBI). Barcelona, Spain, April 5-9, 1998), utilizing phase-sensitive detec- 
tors. With respect to Fig. 1, demodulation is an integral part of the voltmeters 30, 34. Alternatively, demodulation is 
performed by utilizing digital correlation technique, which is accomplished, for example, by the Digital Signal Processor 
40 (Osypka et al. Determination of electrical impedances of tissue at a frequency range of 5 Hz to 20 KHz by digital 
correlation technique. 

[0039] Proceedings of the V. Mediterranean Conference on Medical and Biological Engineering (MEDICON). Patras, 
Greece, August 29-September 1 , 1 989). 

[0040] Voltmeter L 30 also obtains the electrocardiogram (ECG) vector V BC , measured between the left-sided sensing 
electrodes 1 6 and 1 8, which is digitized by an analog/digital converter 56 and fed to the Processing Unit 80. Voltmeter 
Ft obtains the ECG vector V FG , measured between the right-sided sensing electrodes 24 and 26, which is digitized by 
an analog/digital converter 60 and fed to the Processor Unit 80. 

[0041 ] It is understood that more ECG vectors can be obtained by paired combination of sensing electrodes between 
the left hemi-thorax (16, 18) and the right hemi-thorax (24, 26). The measurement of additional ECG vectors requires 
additional voltmeters and analog/digital converters (ADC) connected to the Processing Unit. 

[0042] The Processing Unit automatically, or the operator manually, determines the most appropriate ECG vector, 
or superimposes several ECG vectors to achieve a resulting mean, or reference ECG 62. Alternatively, the outputs of 
several voltmeters are fed into a separate multiplexer. The output of this multiplexer is controlled by the operator, or, 
automatically, by the Digital Signal Processor. Alternatively, the ECG obtained from an external ECG monitor can be 
used as the reference ECG 62. 

[0043] The Processing Unit 80 separates the quasi-constant base impedance, Z 0 , from the time-varying, cardiogenic 
change, AZ(t), or, if the admittance approach is used, the quasi-constant base admittance, Y 0 , from the time-varying, 
cardiogenic change, AY(t). Details of the subsequent processing applied are described below with respect to Figures 
9 and 10. 

[0044] The Processing Unit 80 in Figure 1 is divided into a Digital Signal Processor 40 and a Microprocessor 66. 
Here, the microprocessor 66 establishes the interface between the Digital Signal processor 40 and an operator. The 
functions described as being part of the Digital Signal Processor are not limited to the implementation in the Digital 
Signal Processor exclusively but may be implemented in the Microprocessor, orvice versa. Alternatively, the Processing 
Unit may consist of either a Digital Signal Processor or a Microprocessor. 

[0045] The subject's weight, and other data, is entered via a keyboard 68. Alternatively, data is entered via a touch- 
screen 70, or via a digital interface 72. 

[0046] The stroke volume (SI/) is calculated by the Processor Unit 80 according to the preferred formula 



and wherein the parameters used for this calculation are those which have been input via the keyboard 68, the touch- 
screen 70, the digital interface 72 or those which have been determined in the Processing Unit 80 as set out below 
with respect to figures 6-10. The calculated stroke volume, in conjunction with related cardiovascular parameters, is 
then displayed on a numerical or graphical screen 70. Alternatively, or in addition, it is transmitted via the data interface 
72. 

[0047] The employment of a separate AC source for each tetrapolar electrode array allows measurement of the skin- 
electrode impedance, and, thus, monitoring of the skin-electrode contacts of the current electrodes related to this 




(i) 



wherein 
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electrode array. The apparatus according to Fig. 1 has the capability to monitor skin-electrode contacts of the left and 
right hemithorax individually. Qualitatively, a comparator circuit (not shown) can determine whether or notthe AC source 
is overloaded because, for example, of an infinite load (break of connection, loose electrode). Quantitatively, a voltmeter 
can be connected, or temporally switched, to the outputs of the AC source. The voltmeter measures the voltage across 
5 the skin-electrode interfaces and the thorax. The impedance recognized by the AC source is determined as the ratio 
of the voltage measured and the known AC current applied. 

[0048] Fig. 2 schematically shows an alternative embodiment of the apparatus according to the present invention. 
This embodiment is identical with the apparatus as described with respect to Fig. 1 , but with the exception that only a 
single AC source 28, a single Voltmeter 30 and a single tetrapolar surface electrode array are employed. AC Source 
10 L 28 is connected to the current electrodes 12 and 14, and Voltmeter L 30, connected to the sensing electrodes 16 
and 18. By employing a single electrode array alone (in Fig. 2 the left-sided electrode array), the surface ECG vector 
V BC 54 is available. 

[0049] Fig. 3 schematically shows a further embodiment of the apparatus according to the present invention. This 
embodiment is identical with the apparatus described with respect to Fig. 2, but with the exception that the tetrapolar 
15 electrode array is applied across the thorax. AC Source L 30 is connected to the current electrodes 20 and 14, and 
Voltmeter L 30, connected to the sensing electrodes 24 and 18. By employment of the cross-thoracic electrode array 
alone, the surface ECG vector V FC 54 is available. 

[0050] Fig. 4 schematically shows an alternative embodiment of the apparatus according to the present invention. 
This embodiment is identical with the apparatus described with respect to Fig. 1 , but with the exception that only one 
20 AC source, AC Source L 28, is employed, and interfaced, in parallel, to the current electrodes 12, 14 of the left sided 
and to the current electrodes right sided electrode array, 20, 22. The capability of the apparatus as shown in Fig. 4 to 
monitorskin-electrode contactforthe current electrodes for left and right sides individually, is compromised, with respect 
to Fig. 1. 

[0051] Though the above description was related to the measurement of the impedance (or admittance) of the pa- 
25 tient's thorax utilizing surface electrodes, the apparatus and method according to the invention are not restricted to 
this application. In particular, it is also possible to insert electrodes located on a catheter/probe into the esophagus of 
a patient. 

[0052] Fig. 5 schematically shows a further embodiment of the apparatus according to the present invention. This 
apparatus is identical to the apparatus as described with respect to Fig. 2, but utilizing a tetrapolar electrode array 
30 located on an esophageal catheter / probe. AC Source L 28 is connected to the current electrodes, and Voltmeter L 
30 is connected to the sensing electrodes. The esophageal ECG vector V ES and the esophageal impedance signal V E 
are obtained from the voltage sensing electrodes. The impedance (or admittance) changes as a result of the pulsatile 
vessel flow can be measured. 

[0053] If the apparatus shown in Fig. 5 is used, principally, the stroke volume can be calculated according to formula 
35 (1) given above, wherein coefficients and exponents, and the given implementation of m z rSLS described below have 
to be adapted. The use of electrodes inserted into the esophagus is, for example, disclosed in US patent 4,836,214, 
which is incorporated herein by reference. 

[0054] In the following, it is described how the various parameters used in formula (1) above are obtained. 
[0055] Fig. 6 illustrates the variation of thoracic impedance with a ventilation cycle and with each heart beat. The 
40 electrocardiogram (ECG) on top is a reference for the impedance changes related to the cardiac cycle. The major part 
of thoracic impedance, the base impedance, Z 0 , is obtained as the moving average of measured thoracic impedance 
over a period of, for example, 5 seconds. In a normal healthy subject, Z 0 is approximately 30 £2, and does not change 
from beat to beat. Su-perimposed on Z 0 are changes in impedance (AZ) corresponding to both respiration (AZ RESP ) 
and pulsatile blood flow (AZ CARD ). The total thoracic impedance at any time thus equals: 

45 

Z(t) = Z 0 + AZ(0 = Z 0 + AZ RESP (t)+AZ CARDfAC (t). 

[0056] In Fig. 6, the respiration cycle begins with maximal inspiration, where air in the lungs causes an increase in 
50 thoracic impedance, compared to the base impedance, Z 0 . During expiration, the ratio of air to liquid in the thorax 
decreases, as does the thoracic impedance. 

[0057] If respiration is suppressed, or the corresponding effect on the impedance (or admittance) signal is filtered 
out, only the cardiac-induced pulsatile impedance component, AZ CARDlAC , remains (Fig. 7). For simplicity, in the fol- 
lowing, AZ is referred to the cardiac-induced impedance change, that is, the impedance change due to ventricular 
55 ejection, AZ CARDIAC . 

[0058] Fig. 8 contains parallel tracings of a scalar surface electrocardiogram (ECG) 100, the thoracic cardiogenic 
impedance pulse, AZ(t) 102, and the rate of change of thoracic impedance, dZ(t)/dt 104. 

[0059] The sequential, nearly synchronous, electrical depolarization of the atrial and ventricular heart muscle cham- 
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bers can be electrically sensed and displayed, and the electrical waveform is characterized, by accepted convention, 
as the 'PQRST 1 complex within the ECG tracing. The 'PQRST complex includes the P-wave, corresponding to the 
atrial depolarization wave, the Q-wave (labeled 'Q'), corresponding to depolarization of the inter-ventricular septum, 
the R-wave (labeled 'R'), corresponding to ventricular chamber muscle depolarization, and the T-wave, which repre- 
sents the re-polarization of the ventricular myocardial cells. 

[0060] The R wave is determined, for example, by continuously obtaining the ECG signal amplitude and, by process- 
ing its first time-derivative, its slope. In the event that the absolute slope of the ECG signal exceeds a certain threshold, 
the ECG Processor 62 establishes a time window, within which the detected absolute peak of the ECG signal represents 
the temporal occurrence of the R wave. 

[0061] The time interval between two consecutive R waves is defined as the R-R interval (T RR , Fig. 8). In the appa- 
ratus according to the preferred embodiment, the R-R interval represents the cardiac cycle period. However, alterna- 
tively, other intervals such as, for example, the Q-Q interval can be used to determine the cardiac cycle period within 
the scope of the invention. 

[0062] The R wave peak magnitude of each 'PQRST' complex serves as the temporal reference for the processing 
of AZ(t) and dZ(t)/dt. The point Q precedes the peak R wave by approximately 50 ms and is referred to as the onset 
of electro-mechanical systole. The time interval between point Q and the opening of aortic valve (point labeled 'B') is 
known as the pre-ejection period, T PE . The time interval between point B and the closure of the aortic valve (point 
labeled 'X') is defined as left-ventricular ejection time, T LVE . The point labeled 'C indicates the maximal rate of decrease 
of impedance, i.e. a minimum of dZ(t) /dt.The nadir of dZ(t)/dtai this point in time is further referred to as 



[0063] The point labeled 'Y' represents the temporal occurrence of pulmonic valve closure. The point labeled 'O' 
occurs in diastole and is known to correspond to the early phase of rapid ventricular filling. 

[0064] The value of T LVE can be automatically determined by a computer analysis in which predetermined criteria 
are used. The criterion for the determination of point B is the steep decrease of dZ(t)/dt afterwards, whereas point X 
is the first peak of dZ(t)/dt following 



The determination of these two points of interest is, hence, easy to perform by someone skilled in the computer art. 
The latter automatic method is illustrated in Fig. 9 by the means for determining T LVE 2~\8 as part of the Processing 
Unit (80) and, more specifically, of the Signal Processor (40). Alternatively, the value of T LVE can be determined manually 
by the operator and then, via the keyboard 68, entered into computation. 

[0065] The block diagram of Fig. 9 illustrates the acquisition of the ECG and the thoracic impedance Z(t), with em- 
phasis, in particular, on the Digital Signal Processor (DSP) 40 as part of the Processing Unit 80 (indicated by dashed 
lines). The human subject 10 is shown schematically. The voltage controlled AC source 200 applies an alternating 
current l AC to the subject's thorax. The voltage controlling the current source is generated externally to the Processing 
Unit, or by a synthesizer 202 integrated into the Processing Unit (indicated by dashed line from 202 to 200). In the 
preferred embodiment, the AC source provides a current of constant magnitude, independent of load within reasonable 
limits. Then the AC magnitude must be made known to the Processing Unit 80. In the event that the magnitude of the 
AC source is not constant, it must be measured and recognized by the Processing Unit, as indicated by the dashed 
arrow from 200 to 210. 

[0066] Because of l AC applied to the thorax, a voltmeter 204 can measure the voltage U AC . This voltage contains 
the signal proportional the applied AC and the (unknown) thoracic impedance carrier signal, modulated on a carrier 
frequency, and the ECG signal obtained between the sensing electrodes. Within the voltmeter 204, filters are utilized 
to separate the ECG signal from the applied AC related signal. 

[0067] The ECG signal is the input to an ECG unit 62, which determines the ECG signal used for temporal reference 
for the impedance processing. Alternatively, or in addition, the ECG can be recorded and processed by a source 206 
separate and external to the apparatus described herein (indicated by dashed arrows from 10 via 206 to 62). An 
apparatus 208 determines the RR interval, T RR , from the reference ECG 100 (see also Fig. 6) provided by the ECG 
processing unit 62. Alternatively, T RR can be determined as the time interval between two consecutive occurrences of 
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MIN 



, which approximates T RR (indicated in Fig. 8 and, by a dashed line at 230, in Fig. 9). In an alternative embodiment of 
the invention, the units 62 and 208 are not part of the processing unit 80 but are external devices. 
[0068] The voltmeter 204 eliminates, by demodulation, the AC carrier frequency from the portion of U AC correspond- 
ing to the applied AC. The apparatus 210 determines the impedance Z(t) by calculating the ratio of voltage obtained 
and alternating current l AC applied: 



[0069] A low-pass filter (LPF) 212 is applied to Z(t) in order to obtain the base impedance, Z 0 . A high-pass filter 
(HPF) 214 is applied to Z(t) in order to obtain the thoracic cardiogenic impedance pulse, AZ(t) 102 (see also Fig. 6). 
The purpose of the high-pass filter is also to eliminate impedance changes due to pulmonary ventilation (respiration). 
A differentiator 216 determines the first time-derivative, or slope of AZ(t), that is, dZ(t)/dt. It is also referred to as the 
rate of change of thoracic impedance 104 (see also Fig. 6). 

[0070] With reference to Fig. 6: The dZ(t)idt signal 104 exhibits characteristic landmarks, as described. The left- 
ventricular ejection time T LVE is determined from the dZ(t)/dt signal. Applying basic curve mathematical discussion, 
one skilled in the art can identify the temporal occurrence of aortic valve opening, point B (see arrow), as the "notch" 
just before the steep down-slope of dZ(t)/dt (after the R wave, but prior the point C). Aortic valve closure, labeled as 
point X and pointed to with an arrow, corresponds to a dZ(t)/dt peak after point C. The digital signal processor (DSP) 
40 obtains these points B, C and X, automatically from a processing unit 21 8 (Fig. 7). This unit determines T LVE as the 
time interval between point B and point X. 

[0071] Turning back to Fig. 9, a peak detector 220 is applied to the dZ(t)/dt signal 104 in order to obtain the peak 
rate of change of impedance during systole, see point C in Fig. 7, and its occurrence in time. The ECG provided by 
unit 62 is utilized as a temporal reference. The output of the peak detector, relevant for the SV determination, is the 
absolute peak rate of change of impedance, 



[0072] The left-ventricular ejection time, T LVE , the RR interval, T RR , the base impedance, Z 0 , and the absolute peak 
rate of change of impedance, 



, are furthermore transferred to the microprocessor (uP) 66. The uP 66 determines from the parameters, measured 
and processed by the DSP 40, and other parameters entered, for example, via the keyboard 68 (specifically: weight) 
the stroke volume (SV). The display 70 connected with the uP 66 illustrates the SV and the values of other related 
cardiodynamic parameters. Alternatively, a touch screen can be implemented instead of a display, enabling the operator 
to enter weight and other demographic data via the screen. The iiP 66 can receive data obtained by other, external 
devices, for example, T RR and/or T LVE , through a data interface 72, or send data to other, external devices, such as 
patient monitors. 

[0073] The block diagram of Fig. 10 illustrates the acquisition of the ECG and the thoracic impedance Z(t), with 
emphasis, in particular, on the microprocessor (uP) 66. The uP 66 receives the value for Z 0 from the DSP 40. A unit 
300 calculates the reciprocal of Z 0 , which is then multiplied 302 with the value of 




dzjty 

dt j 
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dm 

dt 



received from the DSP. This product is applied to unit 304, which determines the square root. The result equals 




[0074] which is an integral input for the SV calculation. 

[0075] The uP utilizes the value for Z 0 to determine the index of transthoracic specific impedance 308, further referred 
to as m Zf . This index reflects the presence or absence of abnormal lung water, and is within the scope of the invention. 
m Zf relates to the magnitude or degree of abnormal shunting or bypassing of applied AC around the V EFF , via addi- 
tional abnormal conductive pathways. The critical level of base impedance is defined as Z c , where Z c is greater than 
15 Q. and less than 25 £2, i.e. 15 Q, < Z c < 25 O. In the preferred embodiment, Z c = 20 n (Critchley LAH et al. The 
effect of lung injury and excessive lung fluid on impedance cardiac output measurements in the critically ill. Intensive 
Care Med 2000; 26: 679-685; Critchley LAH et al. Lung fluid and impedance cardiography. Anesthesia 1998; 53: 
369-372; Shoemaker WC et al. Multicenter study of noninvasive systems as alternatives to invasive monitoring of 
acutely ill emergency patients. Chest 1998; 114: 1643-1652; Shoemaker WC et al. Multicenter trial of a new thoracic 
electrical bioimpedance device for cardiac output estimation. Crit Care Med 1994; 22: 1907-1912). 
[0076] In the normal cardiopulmonary state, indicated by Z 0 > Z c , mZ f equals 1. In the abnormal cardiopulmonary 
state, i.e. in the presence of excess thoracic liquids (Z 0 <Z C ), mZ f is less than 1 and greater than 0, i.e. 0 < mZ f < 
1 , and is calculated accordingly: 

mZ = jzS-z c z 0 +c 2 



2Z C 2 +Z 0 2 -3Z C Z 0 + C 2 



[0077] C 2 is a constant and, in the preferred embodiment, taken to 0. In a simplified version of the invention, m Z f is 
taken to be 1 for all values of Z 0 . 

[0078] A unit 310 calculates the reciprocal value of mZ f t or a power of it. The output of 310 is multiplied in unit 312 
with the output of a unit 314 that calculates the mass-based volumetric equivalent of thoracic blood volume in the 
stable, normal cardiopulmonary state. Unit 314 requires the input of the weight of the subject 10 under investigation 
(indicated by the dashed line at 316). In the preferred embodiment, weight is entered via the keyboard 68 or the touch 
screen 70. Alternatively, the value for weight is entered elsewhere and received via a data interface 72 (indicated by 
dashed line at 318). The output of the multiplier 312 is V EFF with 

C -W x 
V = 3 , 

EFF ry N ' 



where C 3 is taken to be 13, but can alternatively have any other value in a range of 0,01-15. In the case of the embod- 
iments shown in Figures 1 -4, C 3 is preferably comprised in the range of 11-15. In the case of the embodiment shown 
in Fig. 5, C 3 is preferably comprised in the range of 0.01 -2.00. In the preferred embodiment, the exponent for weight, 
X, is taken to be 1 .025. with its limits otherwise being 0.9-1 .1 , which are extrapolated from data presented in the article 
by Holt et al. with the title "Ventricular volumes and body weight in mammals", Am. J. Physiol. 1968; 215:704-715. In 
the preferred embodiment, the exponent for mZ f t N, is taken to be 1 .5, with its limits otherwise being 1 .0-2.0. 
[0079] V EFF \s, according to the model used here, the mass-based volumetric equivalent of the thoracic blood volume 
in the stable, normal state. V EFF also represents the total thoracic liquids in unstable cardiopulmonary disease states. 
These conditions are characterized by the abnormal presence of excess thoracic liquids. In the articles by Critchley et 
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al. (Lung fluid and impedance cardiography. Anesthesia 1998; 53: 369-372; The effect of lung injury and excessive 
lung fluid on impedance cardiac output measurements in the critically ill. Intensive Care Med 2000; 26: 679-685) and 
Shoemaker et al. (Multicenter trial of a new thoracic electrical bioimpedance device for cardiac output estimation. Crit 
Care Med 1994; 22: 1907-1912; Multicenter study of noninvasive systems as alternatives to invasive monitoring of 
acutely ill emergency patients. Chest 1 998; 114:1 643-1 652) the impact of excess thoracic liquids related to SV deter- 
mination by means of electrical bioimpedance have been observed. 
[0080] The volume V EFF , determined as 



is an integral part of the preferred embodiment and the new SV equation proposed within. With proper scaling, other 
volumes such as, for example, the ones defined by Sramek and Bernstein, based on weight deviation from ideal weight 
and height, can be used instead. Employment of other volumes is at the expense of accuracy over a wide spectrum 
of subjects, because of their body habita and disease states. 

[0081] The DSP 40 provides the measured values for the left-ventricular ejection time, T LVE , and the RR interval, 
T RR . Alternatively, T LVE can be manually entered via the keyboard 68 (indicated by dashed line at 320), or measured, 
or entered elsewhere and transmitted via a data interface 72 (indicated by dashed line at 322). 
[0082] Alternatively, T RR can be manually entered via a keyboard 68 (indicated by dashed line at 324), or measured, 
or entered elsewhere and transmitted via a data interface 72 (indicated by dashed line at 326). Unit 328 determines 
the reciprocal value of T RR} which equals the human circulatory system frequency: 




[0083] The circulatory frequency, f 0 , or its reciprocal, T RR , can be averaged for a plurality of periods. For example, 
these values can be averaged over the previous ten cardiac cycles ("moving average"). 

[0084] Alternatively, the circulatory frequency, f 0 . can be entered manually by the operator trough the keyboard 68, 
or transmitted via a data interface 72. 

[0085] The value for heart rate (HR, in beats per minute) is calculated by multiplying 334 the circulatory system f 0 
with 60. Alternatively, the heart rate, HR, can be entered manually by the operatortrough the keyboard 68, or transmitted 
via a data interface 72. 

[0086] Unit 330 determines the square root. The output of unit 330 is the Bazett transformation (Bazett MC. An 
analysis of the time relations of electrocardiograms. Heart 1920, 7: 353-364), which, when multiplied in unit 332 with 
T LVE , normalizes T LVE for system mechanical frequency. This normalized T LVE , also known as corrected flow time, FT C , 



is an integral part of the preferred embodiment and the new SV equation proposed therewith . Other embodiments may 
use T[_\/ E instead of FT C for the SKcalculation, at the expense of accuracy at higher heart rates. 
[0087] A multiplier 336 calculates the product of V EFF 



v, 



liFF ~ 



c 2 -w x 



r 



FT,. = 




/ 




and FT C , which equals the SV approximated by the apparatus ([SV\ = mL): 
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sv=v> 



EFF 




= V, 



EFF 




[0088] The value of SV is displayed on a numerical or graphical screen 70. Alternatively, or in addition, it can also 
be transmitted to a data interface 72. 

[0089] Cardiac output [L/min] is then calculated from Sl/[mL] as 



co = 



SV 60 
T„ 1000 



[0090] The measured values can be averaged over a plurality of periods. For example, these values can be averaged 
over the previous ten cardiac cycles ("moving average"). 

[0091] Although Fig. 1-5 and 7-8 indicate that the majority of the functional units are implemented into a processing 
unit, namely a signal processor and a microprocessor, part or all of the functions can be arranged as individual circuit- 
ries. 

[0092] Furthermore, the approximation of SV according to this invention is not limited to the impedance method, but 
can be performed using the admittance approach. With 



Y(t)= 



1 

z(ty 



and 



dY(ty \ ^± 



dZjf) 

dt 



the SV \s approximated according to 



SV = V E 



EFF 



( dY(t) \ 

I dt J 



MAX t t pf 



With respect to the admittance approach, V EFF is determined by 



V EFF =C 3 -W x .mY T \ 



where C 3 is taken to be 1 3 ; but can alternatively have any other value in a range of 0.01-1 5. In the preferred embodiment, 
the exponent for weight, X, is taken to be 1 .025, with its limits otherwise being 0.9-1 .1 . In the preferred embodiment, 
the exponent for m Z f , N, is taken to be 1 .5, with its limits otherwise being 1 .0-2.0. 

The value for Y 0 is utilized to determine the index of transthoracic specific admittance, or conductivity, further referred 
to as m Z 7 . This index reflects the presence or absence of abnormal lung water, and is within the scope of the invention. 
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m Zf relates to the magnitude or degree of abnormal shunting or bypassing of applied AC around the V EFF , via addi- 
tional abnormal conductive pathways. The critical level of base admittance is defined as Y c , where Y c is greater than 

0. 04 £2~ 1 (corresponding to 25 £2) and less than 0.0667 £2" 1 (corresponding to 15 Q), 

1. e. 0.04 Q" 1 < Y c < 0.0667 a- 1 . In the preferred embodiment, V c =0.05 (corresponding to 20 Q). 

In the normal cardiopulmonary state, indicated by Y 0 < Y c , mZ? equals 1. In the abnormal cardiopulmonary state, i. 
e. in the presence of excess thoracic liquids (Y 0 > Y c ), m Z T is greater than 1 , i.e. m Z T > 1 , and is calculated accord- 
ingly: 



C 2 is a constant and, in the preferred embodiment, taken to 0. In a simplified version of the invention, m Z f is taken 
to be 1 for all values of Y 0 . 

[0093] It is noted that, when electrical admittance is determined, instead of electrical impedance, the processing 
performed by the DSP 40 is similar. In this case, the DSP obtains the base admittance, Y 0 , after applying a low-pass 
filter to Y(t), which is the ratio of l AC 200 to U AC 204: 



[0094] The application of a high-pass filter to AY(t) and a differentiator reveals the rate of change of cardiogenic 
admittance, dY(t)/dt. In fact, the dY(t) /dt signal appears, in approximation, as an inverted dZ(t)/dt waveform. In the 
case of the admittance approach, the peak detector determines the peak rate of change of admittance, 



1. Apparatus for determining an approximate value for a stroke volume SV (in milliliter) of a subject's heart, comprising 

a) means for measuring an electrical impedance Z(t) of a part of the subject's body, wherein a value of said 
electrical impedance Z(t) changes with time t as a consequence of the beating of the heart; 

b) means for determining a base impedance Z 0 as a part of said electrical impedance Z(t) which does not 
change significantly during a period of one cardiac cycle; 

c) means for determining a peak magnitude 



of a temporal derivative — ^ of said electrical impedance Z(t), indicating an absolute maximum rate of change 
of said electrical impedance Z(t) during a systolic period of the cardiac cycle; 

d) means for determining a left ventricular ejection time, T LVE ; 

e) means for determining the cardiac cycle period T RR of the heart; and 

f) means for calculating said approximate value of the stroke volume SV wherein said calculating means is 
adapted to evaluate a formula 



2Y Q 2 +Y C 2 -3Y C Y 0 +C 2 
Y 0 — Y c Y 0 + C 2 




Claims 
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sv = v EFF .cv 



dZ(t) 

dt 



MIN 



•T. 



LVE 



wherein 0.15 < n < 0.8 and 0 < m < 1 .5, and wherein V EFF is an approximate value of the subject's volume of 
electrically participating tissue, and wherein C 1 is a constant including 1 . 

2. Apparatus of claim 1 , wherein 0.3 < n < 0.65 ; especially 0.45 < n < 0.55. 

3. Apparatus of claim 1 , wherein m=1 - n, especially m = n= 0.5. 

4. Apparatus of claim 1 , wherein 

c 3 -w x 

V EFF = ^ (in milliliter), 



wherein W is the subject's weight in kilogram (kg), 

wherein C 3 is a coefficient with constant value and X and N are exponents with constant values, 
wherein m Z T =1 for Z 0 > Z c , and 

mZ= = — z c 2 -z c z 0 +c 2 — 

2Z C 2 +Z 0 2 -3Z C Z 0 +C 2 

for 

z 0 <z c , 

and wherein Z c is a constant, and wherein C 2 is a constant including 0. 

5. Apparatus of claim 1 , wherein V EFF = C 3 -W x (in milliliter), wherein W is the subject's weight in kilogram, and 
wherein C 3 is a coefficient with a constant value and X is an exponent with a constant value. 

6. Apparatus of claim 4 or 5, wherein C 3 is a value in the range of 0.01 - 1 5, especially approximately 13. 

7. Apparatus of claim 4 or 5, wherein X is a value in the range of 0.9 -1.1, especially approximately 1 .025. 

8. Apparatus of claim 4, wherein N is a value in the range of 1 .0 - 2.0, especially approximately 1 .5. 

9. Apparatus of claim 4, wherein Z c is a value in the range of 15 - 25 £1, especially approximately 20 Q. 

10. Apparatus of any of claims 1 to 9, wherein said means for measuring said electrical impedance Z(t) comprises: 



at least two pairs of electrodes; 

a current source generating an alternating current l(t) of predetermined amplitude; 
wherein one pair of electrodes is adapted to be connected to said current source; 
means for measuring a voltage U(t) caused by applying said alternating current; 
wherein one pair of electrodes is connected to said means for measuring the voltage U(t); 
means for determining said electrical impedance Z(t) from the voltage U(t) and the current l(t). 
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1 1 . Apparatus of any of claims 1 to 10, wherein said means for determining the peak magnitude 



f dz(tn 

I dt J 



comprises: 

means for determining AZ(t) from Z(t); 

means for calculating for at least the systolic period of one cardiac cycle; 
means for determining tine maximum of an input function. 

12. Apparatus of claim 1 1 , wherein said means for determining AZ(t) is a high-pass filter. 

13. Apparatus of claim 1 , wherein said means for determining the base impedance Z 0 is a low-pass filter. 

14. Apparatus of any of claims 1 to 13, wherein said means for determining T LVE determines T LVE by determining by 
analysis of 

dt 

a point in time when an aortic valve opens; 

a point in time when the aortic valve closes; and by calculating a time difference of said closing point in time 
and said opening point in time. 

15. Apparatus of any of claims 1 to 14, wherein said means for determining the cardiac cycle period T RR comprises 
means for analyzing at least one of a group of Z(t), AZ(t), and 



f dZjtT j 
I dt J 



16. Apparatus of any of claims 1 to 14, wherein said means for determining the cardiac cycle period T RR comprises 
means for measuring an electrocardiogram and means for analyzing the measured values. 

17. Apparatus of any of claims 1 to 1 6, wherein at least one of said means for determining 



f dZ(t) ' 
I dt J 



Z 0 , T LVE and T RR and said means for calculating are comprised in a processing unit. 

18. Apparatus of any of claims 1 to 1 7, further comprising means for calculating an approximate value for a cardiac 
output CO of the subject's heart (in liter/minute) : wherein said calculating means is adapted to evaluate a formula 



CO = SV- 



1 60 



T RR 1000 



19. An apparatus for determining an approximate value for a stroke volume SV of a subject's heart, comprising 

a) a current source outputting an alternating current l(t) of predetermined amplitude to two electrodes; 

b) means for measuring a voltage U(t) caused by said alternating current between two electrodes; and 

c) a processing unit receiving at least a signal representative of U(t), said processing unit being adapted to: 
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calculate an impedance Z(t) from the voltage U(t) and a value of the current l(t); 
input Z(t) into a low-pass filter, an output of said low-pass filter being Z 0 ; 
input Z(t) into a high-pass filter, an output of said high-pass filter being AZ(t); 
calculate a peak magnitude 



f dZ(t) 
I dt 



of 



dZ(t) , 
dt ' 

determine a left ventricular ejection time T LVE from at least one of Z(t), AZ(t) and by using predeter- 
mined criteria; 

determine a cardiac cycle period T RR of the heart from at least one of Z(t), AZ(t) and by using prede- 
termined criteria; 

calculate SV according to a formula 



sv = v EFF c 1 



dZ(t)> 
dt Jmin 



^ 1 r 



LVE 



RR J 



wherein 0.15 < n < 0.8 and 0 < m < 1 .5, especially n = m = 0.5 : 

and wherein V EFF is an approximate value of the subject's volume of electrically participating tissue, and wherein 
C 1 is a constant including 1 . 

20. An apparatus for determining an approximate value for the stroke volume SV of a subject's heart, comprising 

a) a current source outputting an alternating current l(t) of predetermined amplitude to two electrodes; 

b) means for measuring a voltage U(t) caused by said alternating current between two electrodes; 

c) means for measuring an electrocardiogram; and 

d) a processing unit receiving at least a signal representative of U(t) and measured values of said electrocar- 
diogram, said processing unit being adapted to: 



calculate an impedance Z(t) from the voltage U(t) and a value of the current l(t); 
input Z(t) into a low-pass filter, an output of said low-pass filter being Z 0 ; 
input Z(t) into a high-pass filter, an output of said high-pass filter being AZ(t); 
calculate a peak magnitude 



f dZ(tn 
I dt J, 



MIN 



Of 



dZ(t) . 



determine a left ventricular ejection time T LVE from at least one of Z(t), AZ(t) and by using predeter- 
mined criteria; 

determine a cardiac cycle period T RR of the heart from the measured values of said electrocardiogram; 
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calculate SV according to a formula 



sv = v EFF -c 1 



fdZ(t) 
I dt 



MIN 



r 1 \ m 



vTrr J 



'LVE 



wherein 0.15 < n < 0.8 and 0 < m < 1 .5, especially n=m=0.5, and wherein V EFF is an approximate value of the 
subject's volume of electrically participating tissue, and wherein C-| is a constant including 1 . 

21 . An apparatus for determining an approximate value for a stroke volume SV of a subject's heart, comprising 

a) a current source outputting an alternating current l(t) of predetermined amplitude to two electrodes; 

b) means for measuring a voltage U(t) caused by said alternating current between two electrodes; 

c) means for measuring an electrocardiogram and means for calculating a cardiac cycle period T RR of the 
heart from measured values of said electrocardiogram; 

d) a processing unit receiving at least a signal representative of U(t) and a signal representative of T RR) said 
processing unit being adapted to: 

calculate an impedance Z(t) from the voltage U(t) and a value of the current l(t); 
input Z(t) into a low-pass filter, an output of said low-pass filter being Z 0 ; 
input Z(t) into a high-pass filter, an output of said high-pass filter being AZ(t); 
calculate a peak magnitude 



dZ(t) 
dt 



MIN 



Of 



dZ(t) . 



determine a left ventricular ejection time T LVE from at least one of Z(t), AZ(t) and by using predeter- 
mined criteria; 

calculate SV according to a formula 



sv = v EFF -c 1 



dZ(t} 
dt 



MIN 



r 1 v" 



LVE' 



wherein 0.15 < n < 0.8 and 0 < m < 1 .5, especially n=m=0.5, and wherein V EFF is an approximate value of the 
subject's volume of electrically participating tissue, and wherein C-| is a constant including 1 . 

22. Apparatus of claim 21 , wherein said means for calculating the cardiac cycle period T RR is adapted to output a 
signal representative of T RR to the processing unit. 

23. Apparatus for determining an approximate value for a stroke volume SV of a subject's heart, comprising 

a) a current source outputting an alternating current l(t) of predetermined amplitude to two electrodes; 

b) means for measuring a voltage U(t) caused by said alternating current between two electrodes; 
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c) a processor unit receiving at least a signal representative of U(t), said processing unit being adapted to: 

calculate an impedance Z(t) from the voltage U(t) and a value of the current l(t); 
input Z(t) into a low-pass filter, an output of said low-pass filter being Z 0 ; 
input Z(t) into a high-pass filter, an output of said high-pass filter being AZ(t); 
calculate a peak magnitude 



dZ(t) 
dt 



MIN 



Of 



dZ(t) . 



determine a left ventricular ejection time T LVE trorr\ at least one of Z(t), AZ(t) and^^p by using predeter- 
mined criteria; 

calculate SV according to a formula 



sv=v EFF -c 1 



dZ(t) 
dt 



MIN 



•T, 



LVE 



wherein 0.15 < n < 0.8, 

and wherein V EFF is an approximate value of the subject's volume of electrically participating tissue, and wherein 
C 1 is a constant. 

24. Method of determining an approximate value for a stroke volume SV of a subject's heart, comprising the steps of: 

a) measuring an electrical impedance Z(t) of a part of the subject's body, wherein a value of said impedance 
Z(t) changes with time t as a consequence of the beating of the heart; 

b) determining a mean impedance Zq; 

c) determining a peak magnitude 



dZ(t) 

I dt . 



MIN 



of a derivative of said impedance Z(t) over the time t by using the measured impedance Z(t) for at least 
a systolic period of one cardiac cycle; 

d) determining a left ventricular ejection time, T LVE ; and 

e) determining a cardiac cycle period T RR of the heart; 

f) calculating an approximate value of the stroke volume according to a formula 



sv = v EFF -c 1 



dz(t)- 

dt Jm\N 



f 1 y 



■T, 



LVE 



RR . 



wherein 0.15 < n < 0.8 and 0 < m < 1 .5, 
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and wherein V EFF is an approximate value of the subject's volume of electrically participating tissue, and wherein 
C 1 is a constant including 1 . 

25. Method of claim 24, wherein 0.3 < n < 0.65, especially 0.45 < n < 0.55. 

26. Method of claim 24, wherein m = 1 - n, especially m = n = 0.5. 

27. Method of claim 24, wherein 

C 3 -W x 

V EFF = ^ (in milliliter), 

mZ- 

wherein W is the subject's weight in kilogram (kg), 

wherein C 3 is a coefficient with constant valije^gc^X + gnd N are exponents with constant values, 
wherein m Z f = 1 for Z 0 > Z c , and m Z f = E — E_2 — ? — for 

2Z c 2+Z 0 2 - 3Z cV C 2 

z 0 <z c , 

and wherein Z c is a constant, and wherein C 2 is a constant including 0. 

28. Method of claim 24, wherein V EFF =C 3 -W X (in milliliter), 

wherein W is the subject's weight in kilogram, and wherein C 3 is a coefficient with a constant value and X is an 
exponent with constant value. 

29. Method of claim 27 or 28, wherein C 3 is a value in the range of 0.01 - 15, especially approximately 13. 

30. Method of claim 27 or 28, wherein X is a value in the range of 0.9 -1.1, especially approximately 1 .025. 

31. Method of claim 27 or 28, wherein N is a value comprised in the range of 1 .0 - 2.0, especially approximately 1 .5. 

32. Method of claim 27, wherein Z c is a value in the range of 15 - 25 Q, especially approximately 20 n. 

33. Method of claim 24, wherein said impedance Z(t) is measured by applying an alternating current l(t) through the 
part of the subject's body, measuring a voltage drop U(t) in the body caused by the application of said alternating 
current, and calculating said impedance Z(t) according to a formula 

k } i(t) 

34. Method of claim 24, wherein said peak magnitude 

f dZ(tn 
I dt J MIN 

is determined by: 

sending a signal representative of Z(t) through a high-pass filter, an output of said filter being taken to be AZ(t) ; 
calculating — for at least the systolic period of one cardiac cycle; 
determining an absolute magnitude of said derivative. 

35. Method of claim 24, wherein Z 0 is determined by sending a signal representative of Z(t) through a low-pass filter 
an output of said filter being Z 0 . 
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36. Method of claim 24, wherein T LVE is determined by 

determining by analysis of 

a point in time when an aortic valve opens; 

a point in time when the aortic valve closes; and 

and by calculating a time difference of said closing point in time and said opening point in time. 

37. Method of claim 24, wherein the cardiac cycle period T RR is determined by analyzing at least one of the group of 
Z(t), AZ(t),^and 

dt 




for at least two consecutive cardiac cycles. 

38. Method of claim 24, wherein the cardiac cycle period T RR is determined by 

measuring an electrocardiogram, and 
analyzing the measured values. 

39. Method of claim 24, further comprising the step of determining an approximate value for a cardiac output CO of 

the subject's heart, wherein CO (in liter/minute) is determined as the product of SV (in milliliter), T — and a constant: 

'rr 

CO=SV-. 1 60 
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